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Summary. — Recent theoretical results on charged particle interaction with pla- 
narly oriented thin bent crystals are reviewed, with the emphasis on dynamics in 
the continuous potential. Influence of boundary conditions on the volume-reflected 
beam profile is discussed. Basic properties of coherent bremsstrahlung in a bent 
crystal are highlighted. 

PACS 61.85.+P - bent crystal; volume reflection; coherent bremsstrahlung. 



1. Introduction 

Passage of charged particles through a bent crystal, even in a planar orientation, 
is a complex phenomenon, theoretical treatment of which used to rely on computer 
simulation. However, sometimes complexity begets simplicity, which may permit analytic 
advances - like it was commonly practiced before the advent of computers. Actually, 
computer and analytic calculations complement one another: analytic formulas provide 
general understanding of the phenomenon, and may serve for experiment planning, while 
computer can offer precise predictions for an established experimental setting. 

An example of a problem where analytic approach did bring fruit is fast charged par- 
ticle passage through a planarly oriented uniformly bent crystal. Even though the planar 
orientation and uniformity of the bending enormously simplify the dynamics, reducing 
it to radial Id, yet there are impediments for solution of the whole beam-crystal inter- 
action problem: how to analytically describe particle passage through many inter-planar 
intervals, and how to analytically average over the initial particle impact parameters. 
Alleviation came from the use of a simplified (parabolic) model for inter-planar contin- 
uous potential, which made analytic solution for single volume reflection (VR) feasible 
[1] . It appears that precise inter-planar potential shape is of minor consequence anyway, 
and the piecewise harmonic potential model yields fair agreement with the experiment. 

A different story is the bremsstrahlung emitted by the fast particle in a bent crystal. 
There, a large contribution comes from the dipole coherent bremsstrahlung, which can 
be described for an arbitrary inter-planar potential. Analytic theory of dipole coherent 
radiation in bent crystals (CBBC) was constructed in [2]. 
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The present short note surveys the physical picture of volume reflection and radiation 
in a bent crystal, in the pure continuous potential. 

2. Volume reflection 

21. Essence of the effect. - VR effect manifests itself in the same region of E/R ratio 
as channeling, with the proviso that the particle entrance angle with respect to active 
crystallographic planes must be well above critical. However, the correlation between the 
particle deflection angle and the crystal curvature for VR appears to be more intricate 
than that for channeling: 

1. in VR particles deflect to the side opposite to that of the crystal bending; 

2. the deflection angle tends to a finite limit with the increase of R/Rc ratio, although 
R = oo corresponds to a straight crystal and no net deflection; this limiting value is 
of the order of critical angle 9c both for positive and for negative particles, though 
with different numerical coefficients. 




Fig. 1. - Left-side frame - trajectory filling of an exemplary inter-planar interval by an initially 
transversely uniform particle beam (top - positive particles, bottom - negative particles). The 
deficit of positively directed force in VR area causes the negative net deflection of the beam. 
Right-side frame - comparison of ultra- high energy (7? <^ 47?c, fepc 2ttt) with moderately 
high energy {R ^ 47ic, 2ttt <g Iefc) particle passage. In the latter case the particle net 
deflection angle is non-zero and ~ 9c, regardless of the particle charge sign. 



The qualitative explanation of the opposite deflection direction is given in [3]: "For 
a small crystal curvature, that is i? ^ Rc, the turning points of all particles gather 
in a narrow region near the inner wall of a planar channel. The strong electric fleld 
of the crystal plane is directed along R and at the turning points it imparts to the 
particle an angular deflection towards the opposite direction with respect to the crystal 
bending, producing volume reflection." For explanation of feature 2 the latter argument 
is not sufficient yet. It is known that for negative particles the force in the reflection 
point is much smaller than for positive ones, but nevertheless, the deflection angle for 
these cases is of the same order. For a closer look, let us draw a graph of the family 
of particle trajectories in an exemplary inter-planar interval (Fig. 1). The figure reveals 
that in inter-planar channels there emerge regions devoid of particles. This entails deficit 
of force of a definite (positive) sign acting on the beam as a whole; hence, the mean 
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deflection angle must be negative. At sufficiently large R, the void transverse dimension 
is ~ d, so F ~ i^max, and the longitudinal dimension is r = ^/R^A/2 (the channeling 
oscillation timescale)("'^); none of these parameters involve R, therefore, there is a limiting 
value limfl_j.oo {0)^ independent of R: {9)^^ ^ -^F^axT = &c- 

Of course, if one overbends the crystal, all the intra-crystal space becomes uniformly 
covered by the particle flow, then voids and therewith the VR effect disappear. In [1, 2] 
it was proved that ii E oo, i.e., R/Rc 0, then J dhO{h) = 0. 

Ex adverse the existence of VR effect can be explained more concisely: When a high- 
energy particle passes through a bent crystal, after its entrance the angle between its 
velocity and the active planes decreases. When this angle becomes ~ 9c, the particle 
experiences strong influence of the continuous potential fleld on the length r. If this 
length is shorter than ^efc = V'2Rd, it would be capable of making the particle channeled. 
But for motion in a stationary potential, an unbound particle cannot pass into bound 

2 

state. Hence, at ttt < Ibfc, i-C- R > ^Rc, the particle net inclination angle can 
not become smaller than 9c, which implies that the particle will reflect when the angle 
becomes ^ 9c - see right-side frame of Fig. 1. 

2'2. Boundary condition sensitivity . - Although the term 'volume' reflection signifles 
that the deflection arises somewhere in the crystal volume, and therefore should not 
depend on the crystal boundaries, but to some degree, the boundary effects still manifest 
themselves when investigating the final beam profile detail. 




a) b) c) 



Fig. 2. - Possible boundary conditions for VR in a uniformly bent crystal. Left-side frame, 
cases a) and b) (the figure is taken from [4]) - particle entrance through a flat front face of a 
thin bent crystal. At modern practice, the a;-dimension of the crystal is actually much wider 
than its ^-dimension. Right-side frame (figure taken from [5]) - particle entrance through the 
curved lateral surface of a longitudinally extended crystal; this type of boundary conditions is 
also adopted in [6]. 

First of all, even for a uniformly bent crystal, there are several options for particle 
beam entrance to the crystal - see Fig. 2. At practice one usually deals with case b), when 
the particle entrance angle 9q is much greater than critical; the corresponding theory 
was worked out in [7, 1]. But at channeling experiments, there also arises case a), which 



(^) Note that this scale is R/Rc times shorter than scale R9c within which the particle trajectory 
substantially differs from a straight line. At 7? ^ Rc, in principle one should mind the existence 
of two longitudinal scales in the volume reflection phenomenon. 
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was studied in the pioneering work [4]. In the latter case, the particles enter the crystal 
tangentially to the active planes, so, at JS ^ i?c most of the particles are channeled, but a 
small fraction of particles hitting one-sided vicinity of the top of potential barrier belongs 
to volume reflection. For the latter fraction, it is supposed that owing to symmetry of 
the trajectory in a central field with respect to the reflection (minimal-radius) point, 
the deflection angle 9 for case a) is half the value for case b). But yet there is a slight 
dependence of the deflection angle on the impact parameter, and it is this dependence 
which determines the final beam shape. So, a question remains, whether the final beam 
profiles for cases a) and b) are similar. 



R/27rR,. . 



e^<dw/de>, 



00 -dependent 
rainbow spike 




2 2 K 

Fig. 3. - Left-side frame - angular distribution in a VR beam for positive particles, at 7? > 
4i?c, without account of multiple scattering. Solid curve ~ distribution averaged over a small 
variations of initial large angle 6o, dashed curve - for a definite do- Right-side frame - the 
angular distribution for negative particles, plotted for R — 25Rc- Black curve - for boundary 
condition of type b) (averaged over a small vicinity of large angle 6o), gray curve - for case a) 
and stretched in 9 by the factor of 2. The peaks on the figures are due to rainbow phenomena[l] . 



To establish correspondence between the profiles for cases a) and b), neglecting mean- 
while the multiple scattering, note that 9 depends on the impact parameter b only through 
the transverse energy E±_ + Vce{-b), with Vcsi-b) = V{-b) + whereas 9{E±_) 

is a universal periodic function of E±_ with the period AE±_ = ~ "^K) ^ Vq [7]. 
However, dependencies E±(b) are different under different boundary conditions. For 
case b), 9{E±) spans many periods as b varies from — | to |, and E±{b) dependence, 
basically, is linear within each period of 9{E±). But for case a) 9{E±) spans only one 
period, and the dependence of Veg near its top is essential. Now, for positive particles 
this behavior is close to linear, too, and so for positive particles boundary conditions a) 
and b) yield basically isomorphic profiles. On the contrary, for negative particles Vcs{b) 
behavior near the top is quadratic. Thus, the relation between the profiles appears to 
be 20°) (y^) = 9^\u)^ where v is the transverse kinetic energy on the entrance to the 
reflection inter-planar interval, rescaled to vary from to 1. Sec Fig. 3, right-side frame. 

As for boundary condition of type c) (see Fig. 2), which simplifies the theoretical 
problem by making it entirely centrally-symmetric, it yields^ similar to the case b), 
though indicatrix 9{b) has a somewhat different double-periodic behavior [5]. 



(^) In paper [I], footnote 14, it is asserted that for case b) there is no rainbow, but more 
precisely, the rainbow significantly attenuates, though in principle remains. 
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2'3. Multiple scattering . - Intrinsic VR angular distributions shown in Fig. 3 have an 
interesting edgy structure, for negative particle case being yet asymmetric and containing 
a significant rainbow. But at practice, those distinction features will be obscured by 
multiple scattering. However, the extent of the VR area 27rr = tt\/2TR^ ^ 50^m is 
much shorter than the crystal thickness L ^ 1mm. Thereat, the multiple scattering 
eff'ects accumulate mainly away from the VR point (upstream and downstream of it). 
This allows one to view the aggregate distribution as a convolution of successive scattering 
probabilities [7]: 



where dwv.r./d0 is calculated with the neglect of multiple scattering, and a^^^ is the 
mean square angle of multiple scattering in an amorphous target of the same thickness. 
Consideration beyond approximation (1) will be given elsewhere. 

In principle, Eq. (1) can be inverted; it is noteworthy that Fourier transforms of 
dwreni/dO and dwv.r./dO are proportional, and should both have negative regions. But 
unfortunately, under conditions of binned measurement of dw^cai/dO, only first few mo- 
ments of the angular distribution can be accessed reliably. 

2'4. Moments of the final angular distribution. - The complete analytic calculation 
of dwv.i./dO was made within the model of [1]. According to this solution, the limiting 
value for 9^,-[J6c at R/ARc — > oo amounts 7r/2 for positively charged particles, and 1 
for negatively charged particles. More precisely, this is the outer edge of the angular 
distribution. The deflection angle mean value involves an 0{Rc/R) correction: 
(2) 

^dwi-cai f ,nndwv-r. _ j ^f^c(l-^) for pos. charged particles 

d9 I — ^?c ( 1 — 1-3-^ In j for neg. charged particles 



{9) = J cl99^^ = J d9e- 



(the coefiicient 1.3 for negative particles needs more precise evaluation). For protons, 
prediction (2) was compared with experiment in [1], see also Fig. 4 below. 

The second moment, mean square deflection, under assumption (1) possesses the 
property of additivity, even for non-gaussian dwv.r./d9: 

(3) ^Li « / dei9- {e)f / da ' = aL + -al... 

J-00 J-00 V27rCTam 

where ,. ~ J°°^ d9{9 — (0))^ '^^Je' • positively charged particles, acquiring rectan- 
gular intrinsic final angular distribution at < 1, evaluation of (Tv.r. gives 

TT Rc 

(4) tTv r = ~F~;7^c (for positively charged particles) . 

V3 



^v.r. 



Experiment [3] used procedure (3), i.e., formula fiv.r. = \/ Cj^^ai ^ '^am 1 to determine 
Comparison of experimental data with Eq. (4) is shown in Fig. 4 by dashed line. 

For negatively charged particles, where the final profile is skew, its third moment is 
also non-zero: J d9 (9 ^ {9)f ^ ^ d9 {9 - {e)f ^ 0. 
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Fig. 4. - Deflection angle mean value (Eq. (2), solid line) and rms deviation (Eq. (4), dashed 
line) for a parabolic model of inter-planar potential compared with CERN SPS experimental 
data [3] (protons, £=400 GeV). 



2'5. Conditions for quality VR deflection. - Joint consideration of beam spread in 
the continuous potential, multiple scattering in the silicon material, and the initial beam 
angular divergence has led [1] to a system of conditions needed for practical application 
of VR to beam steering at accelerators: 

L E R /20Afrad 

^ ' 1mm ^ 100 GeV ^ m' V 

where ctq is the initial beam rms divergence. Conditions (5) are fairly well met in exper- 
iments [3]. 

2'6. VR at 6o ~ 9c- - The theoretical description of VR in [7, 1] focussed on the case 
^'o 3> 0c; when the process may be regarded as independent of the crystal boundaries. But 
in experiments, VR is often investigated in parallel with channeling, with a continuous 
passage from one regime to another. In the transition regime, the final beam angular 
distribution must acquire substantial ^o-dependence. Investigation of this dependence is 
important as well (see [8]). 

3. Radiation in thin uniformly bent crystals 

Observation of inelastic processes evoked by a fast particle traveling a bent crystal 
may be used for extraction of information about particle trajectories at real conditions. 
Among possible inelastic processes, the coherent radiation in the near-forward direc- 
tion("^) benefits from existence of a local correspondence between the photon ui and the 
longitudinal coordinate z it is emitted from. This relation is particularly simple (ge- 
ometric rather than dynamic) at large oj, corresponding to z away from the VR area. 
There, the particle deflection from the straight line is only perturbative, opposite in sign 
for electrons and positrons. Here we will confine ourselves to discussion of the latter 
simplest case - coherent bremsstrahlung in a bent crystal (CBBC). 




(^) The direct radiation from the fast passing particle is only significant for electrons and 
positrons, whereas for beam deflection applications the primary interest is on protons. However, 
ultra-relativistic positrons have essentially the same trajectory as protons of equal energy, so 
one indirectly can test proton passage through experiments with positrons. 
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3'1. CBBC intensity as a sum of infinitesimal straight-crystal contributions. - The 
first, analytic, calculation of coherent bremsstrahlung in a bent crystal was undertaken 
in [9] (in appendix of a paper primarily dedicated to channeling radiation in a uniformly 
bent crystal). Fourier decomposition of the bent crystal continuous potential yielded 
Fresnel functions, which describe both volume effects (step- like asymptotics), and the 
oscillations due to edge effects. But at practice usually the edge effects are negligible, 
and the Fresnel functions may be replaced by a Heavyside unit-step function, which led 
to Eq. (A. 16) of [9]: 



1 



A4,{R,L) 



A'lp[R,L) - Ornin Je,„i„ 

with da^ , da^ the radiation differential cross-sections in a straight and in a bent crystal. 
Alp the angle between the particle velocity and active crystallographic planes on the 
crystal edge, and ^min — ^^2^ This representation in form of unweighted 0^- 

averaging is only valid at i? = const. Eq. (6) also assumes symmetric orientation of the 
beam with respect to the particle trajectory (when Aip at the entrance from the crystal 
equals Atjj at the exit), though generalization to an asymmetric case is straightforward. 
In fact, at practice such a generalization may be necessary even if the crystal is oriented 
symmetrically with respect to the beam axis, since non-negligible beam divergence com- 
pared to the crystal half-bending-angle makes the orientation asymmetric for individual 
particles [11]. 

A generic representation (valid at variable R and arbitrary orientation), which can 
be obtained based on the stationary phase approximation [2] (^), reads 



J-L/2 dzduj 



,=ie(2)-«oi 



where ^(z) is the transverse coordinate of any of the bent atomic planes. If one utilizes 
here the known formula for the spectrum of coherent bremsstrahlung in a straight crystal 
(8) 



dzdw 27r4m26l2 £;2"^"""^^V'" 2TT\eJl n'^+^'^ V mrWJ ' 2nVdl ' 2EE 



with gmin('^) — ; E' ^ E — oj, 9x{z) = Ciz) — do, equation (7) turns to explicit 
Eq. (30) of paper [2]. To see the correspondence of Eqs. (6) and (7), one observes that 
in a uniformly bent crystal S,{z) = -^ji + const, so integration variables 9x = \z/R — 9o\ 
and z are linearly related. In the ratio of da^' oc L and Aip cx ^ the crystal thickness 
L cancels out, so ^ oc i?. From the viewpoint of Eqs. (7, 8), the proportionality of the 
radiation spectrum to the crystal bending radius arises as 

/n\ dEcBBC dz 2 

(9) __oc-^M^Zepc; 



(*) The replacement of Fresnel functions by step functions in [9] is equivalent to application of 
the stationary phase approximation in the original Fourier integrals. 
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see also Eq. (14) below. 

For a uniformly bent crystal, further on, it is straightforward to accomplish integration 
over Ox and arrive at Eq. (31) of [2], presenting the spectrum through a cubic polynomial 
function D. Averaging over angular distribution in the initial beam and account of 
multiple scattering in the crystal volume is also feasible analytically. 

3'2. Coherence lengths and typical photon energies. - Evaluation of the spectrum of 
bremsstrahlung produced under a definite force F[t) action on a charged particle often 
proceeds in two steps: (i) Fourier-decomposition of the external force, i.e., distinguishing 
virtual photons absorbed by the particle; (ii) conversion of the force spectrum to the 
radiation spectrum. An important physical scale characterizing each stage is its coherence 
length. In general, real photon emission (Zform) and virtual photon absorption (Zefc) 
coherence lengths need not be related. Indeed, in simplest cases one has ^efc — > 
(stochastic scattering in an amorphous medium), or ^efc ^ (uniform external field), 
but /form = always stays finite, providing for those cases the only coherence length to 
refer to. From the experience of the mentioned simplest cases, it had become common to 
call Zform simply 'coherence length', without specification of the nature of the coherence. 

On the other hand, if the bremsstrahlung problem may be treated in the dipole 
approximation, the photon emission coherence is strictly related with the virtual photon 
absorption coherence: it reduces to Lorentz-rescaling of the frequency and integration of 
the intensity over typical ^ 7"^ emission angles (cf. Eq. (11) below). Thence, it suffices 
to analyze the coherence in the external field Fourier decomposition. 

In crystals, ^efc acquires a finite size. In a straight crystal, this is just the distance 
between the planes crossed along the particle path; the irradiation resonance condition 
requires /form to have the same value. Then, 

(10) /form ^ ^EFC - ^76*0 (straight Crystal) . 
At that, the typical radiation energy 

(11) Wo - 27^/-^^ - 7^6*0/^ (straight crystal) 

is just the Lorcntz-rescaled plane-crossing frequency. The radiation spectral intensity is 

estimated as (cf. (8)) 

(12) 

dE _ dE (^^X ^ 

duJ dzduj \^ J ^0 ^form(ti^) 



— = -JTJ-.^ ( — ) -^7 TT' ^form = 9mm(w) < ^'o/o? (straight Crystal) . 



Less trivial situation emerges in a bent crystal. There, the local frequency of plane 
crossing varies along the crystal, whereby there are 2 spatial scales: 

(13) Zefc = V2Rd (bent crystal) , 

on which the integrals from oscillatory gaussians (J (ize*^i?3 ...) converge, and the longi- 
tudinal geometrical scale of the crystal (say, its thickness L) determining the limiting 
frequency of plane crossing. If a tangency point of the particle trajectory with the family 
of bent planes occurs within the crystal, then 

dE f cFV 

(14) — oc — 'ffp (bent crystal). 
du) \ m I 
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Comparing Eqs. (14) and (12), one sees that in both cases the radiation intensity is 
proportional to Z|pQ, which is in line with the understanding [10] of coherence length as 
the length at which radiation amplitudes add up. Besides that, (12) involves a factor 
arising due to spectral overlap of radiation generated within the crystal thickness 
in different coherence intervals (pile-up factor). In a bent crystal, where radiation from 
different coherence intervals does not overlap, factor enters the expression for the 
spectrum extent 

(15) Wend (X min {47r7^L/;|pc, E] 

instead of intensity (14). 

Since /efc ^form appear to be important physical quantities, it is worth estimating 
their typical values. For i? ^ 20 m, d w 2 A, Eq. (13) gives ^efc ^ 100 /zm, whereas 
ratio ^ is usually small (see Eq. (16) below). Therefore, without breaking the 
picture of CBBC, the crystal thickness may be decreased down to fractions of millimeter. 
Lastly, concerning the multiple scattering influence on the radiation, it is quantified by 
parameter |eec. ^ which is larger than LPM parameter , and therefore is more critical. 

3'3. Locality of CBBC generation vs. straight crystal limit. - The observation that 
representation (6) contains a straight crystal limit may appear surprising from the view- 
point that local CBBC theory is based on the stationary phase approximation, whose 
condition assumes 

(16) Zefc < (CBBC locality condition). 

As R, and therewith ^efC: increases, condition (16) must break down. Nevertheless, this 
does not destroy the convergence of the integral, provided 

(17) 4d/|0o| (condition of many- interval crossing in a staight crystal) , 

i. e., the particle crosses a large number of inter-planar intervals even in a straight crystal. 
If (17) holds, with the increase of R one sooner arrives at a condition 

(18) |6'o| > L/2R 

than at (16). That implies that variation of the local plane-crossing frequency in the bent 
crystal is smaller than the frequency mean value, which is equivalent to near straight- 
ness of the crystal, regardless of whether condition (16) holds or not. In [2] it was 
demonstrated that under condition (18), Eqs. (7-8) indeed turn to the familiar formula 
of CBBC in a straight crystal, with the crystal bending radius dropping out. Hence, 
further increase of the radius will be inconsequential. The Zmear-exponent oscillatory 
integral convergence will be achieved at length (10) within the crystal thickness without 
the aid of crystal curvature. Together, conditions (16-17) may be cast into a universal 
expression for the convergence length: Iconv = min {Zefc, 4(i/|0o|} • 

3'4. Radiation polarization. - For planar orientation of the crystal, when the particle 
is mainly subject to the force orthogonal to the active planes, there must be a significant 
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net polarization of the cone of emitted bremsstrahlung, directed orthogonally to the 
planes. According to [?], the polarization degree equals 



where the azimuthal anisotropy parameter N for planar orientation must be close to 1. 
Thereat, the maximal value of Pnet is 50% at w <C -E, while at oj ^ E, E' ^ 0, Pnot — > 0. 

4. Conclusions 

The problem of particle interaction with bent crystals in planar orientation contains 
many opportunities for analytic description. Some successes have already been achieved, 
linear laws of Fig. 4 being an example. Ultimately, for this class of problems there 
is a good perspective to calculate all the relevant observables without MC simulation, 
although the number of problems to deal with is large, and the goal remains remote 
presently. 
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